Microbiome Data Analysis Snakemake Workflows

Integrated Solutions for Better Results

Quick glimpse

• We use the snakemake workflow engine¹ to maintain reproducibility in technical validation, regeneration of results, and improvement of the microbiome bioinformatics analysis.
• We review existing workflows to gain a better insights and regenerate contiguous but related rules.
• Each major step form a separate executable snakemake rule.

We envision to keep fostering on continuous integration and development of highly reproducible workflows.

---

IMAP-Mothur Snakemake workflow

• A mothur-based workflow for microbial profiling using 16S rRNA gene markers.
• Each major step forms a separate executable snakemake rule.
• GitHub repo: https://github.com/tmbuza/imap-mothur

Workflow framework

Microbial profiling classification options

1. Operational Taxonomic Units (OTUs)

• OTUs are clusters of similar sequences and are commonly accepted as analytical units in microbial profiling when using 16S rRNA gene markers.

2. Phylotype

• A phylotype in microbiome research is a DNA sequence or group of sequences sharing more than an arbitrarily chosen level of similarity of a 16S rRNA gene marker.

3. Amplicon Sequence Variant (ASV)

• An ASV in microbiome research is any inferred single DNA sequences recovered from a bioinformatics analysis of 16S rRNA marker genes.

4. Microbial Phylogenies

• Microbial phylogenies are from gene sequence homologies. Models of mutation determine the most-likely evolutionary histories.

---

Appendix

Reference Databases

1. Mothur-based SILVA reference files².
2. Mothur-based RDP reference files³. Note: The RDP database is to classify 16S rRNA gene sequences to the genus level.
3. ZymoBIOMICS Microbial Community Standard (Cat # D6306)⁴. The ZymoBIOMICS Microbial Community DNA Standard is designed to assess bias, errors and other artifacts after the step of nucleic acid purification.

Sample location (demo)

Troubleshooting (in progress)

1. Mothur dist.seqs taking too long.
• Merged reads are too long, probably over 300pb.
• Reads not overlaping when merging the paired reads.
• Too many uniques representative sequences probably caused by lack of overlapping.
• No enough computer power which suggest a use of HPC or Cluster.

References

1. Köster, J., Mölder, F., Jablonski, K. P., Letcher, B., Hall, M. B., Tomkins-Tinch, C. H., Sochat, V., Forster, J., Lee, S., Twardziok, S. O., Kanitz, A., Wilm, A., Holtgrewe, M., Rahmann, S., & Nahnsen, S. (2021). Sustainable data analysis with snakemake. F1000Research, 10. https://doi.org/10.12688/f1000research.29032.2

2. Mothur-based silva reference files. https://mothur.org/wiki/silva_reference_files/

3. Mothur-based RDP reference files. https://mothur.org/wiki/rdp_reference_files/

4. ZymoBIOMICS microbial community DNA standard (cat # D6306). https://www.zymoresearch.com/zymobiomics-community-standard

---

Last updated on January 25, 2023.

These workflows are continuously reviewed, integrated, tested and compiled by Teresia Mrema-Buza.
Github_Repo: https://github.com/tmbuza/imap-mothur.
Github_Page (this document): https://tmbuza.github.io/imap-mothur/.
Visit the Complex Data Insights (CDI) website for more practical user guides (…in progress).